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<3 Abstract 

<N ■ 

^ ■ The mechanism of the exponential transient statistics of Poincare recur- 

rences in the presence of chaos border with its critical structure is studied us- 

' ing two simple models: separatrix map and the kicked rotator ('microtron'). 

ON ■ For the exponential transient to exist the two conditions have been shown 

1 to be crucial: fast (ballistic) relaxation, and a small measure of the critical 
Q | structure. The latter was found to include a new peripheral part (halo) of a 

■ surprisingly large size. First preliminary empirical evidence is presented for a 

. new regime of Poincare recurrences including the transition from exponential 
to exponential statistics. 

1 Introduction: exponential vs. power— law PR 

m ! 

q ■ As is well known any trajectory of a bounded in phase space motion of Hamiltonian 

<§; 

both regular (with discrete spectrum) as well as chaotic (with continuous spectrum) 

r— h . motion. These Poincare recurrences (PR) do not imply a quasiperiodic motion 

C ■ 

which is still a widespread delusion (see, e.g., |J). The difference between regular 

'O : 

and chaotic motions lies in the statistics of recurrences which is usually described 
by the integral distribution P(r) that is by the probability for a recurrence time to 
be larger than r. In a regular motion such a survival probability P(r) has a strict 
upper bound in r while for a chaotic motion r can be arbitrarily long. In both cases 
PR characterize some fluctuations including arbitrarily large ones in chaotic motion. 
The PR statistics proved to be a very powerful and reliable method in the studies 
of chaotic dynamics due to its statistical stability. 



1 Email: chirikov@inp.nsk.su 



1 



To my knowledge, such a method was first used (implicitly) in Ref. |2| for the 
study of a narrow chaotic layer along the separatrix of a nonlinear resonance. The 
result (r > 1) 

P(r) « ^ (1) 

was a surprise as it contradicted the bounded motion in chaotic layer. Indeed, the 
total sojourn time r-P(r) of a trajectory, which is prportional to the measure of the 
chaotic component of the motion, diverges asr-> oo. Later ||, this apparent con- 
tradiction has been resolved simply by increasing r which showed that the exponent 
of the power-law decay also increased from the initial v — 1/2 to v ~ 3/2. 

It is instructive to mention that the origin of a short-time computation in Ref. || 
was in apparently reasonable decision to avoid any rounding-off errors by enormous 
increase of the computation accuracy. As a result, the computation speed, and 
the available motion time, dropped by several orders of magnitude. Generally, for 
exponentially unstable (chaotic) motion such an approach is prohibited whatever 
the computer power. Fortunately, it is also unnecessary for calculating statistical 
characteristics of the motion like P[j) since most of the latter are robust. True, 
the corresponding Anosov theorem was (and can be so far) proved for the very 
simple Anosov systems only. Moreover, such a theorem is even wrong for discontinu- 
ous (discrete) perturbations like rounding-off ones (see, e.g., Refs.||). Nevertheless, 
all the numerical experience confirms a sort of robustness of the statistical behav- 
ior of chaotic systems, at least with some minimal precautions (see, e.g., Refs.0 
for discussion). Notice that without such an 'empirical' robustness the numerical 
experiments with always approximate models would lose any physical meaning! 

A power-law decay -P(t) ~ t~ u , whatever the exponent z/, found in 0, || for a 
bounded motion, was at variance with the exponential decay believed to be a generic 
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case. In Ref . the former was interpreted as a characteristic of a qualitatively new 
structure of the motion near the chaos border in phase space. Later, it was termed 
the critical structure, which was described by a renormalization group |7| (see also 
review || and references therein). 

Since then, the exponential decay has been considered as a property of ergodic 
chaotic motion without any chaos borders. However, in recent numerical exper- 
iments || with an asteroid motion a fairly long transient exponential decay was 
found. Moreover, it persists in the separatrix map also used, just the same map 
which seemed to have been well studied in many previous works (see also || and 
references therein). 

The main purpose of this paper is to reconsider various regimes of PR, and to 
formulate the conditions for their realization using two relatively simple models: 
separatrix and standard maps. Only bounded motion will be considered, with or 
without chaos borders. First, a classical problem of PR in an ergodic system will be 
discussed in some details in Section 2. Then, in Section 3, the analysis of various PR 
regimes in the separatrix map will be presented aimed to resolution of the apparent 
contradiction mentioned above. In Section 4, PR in the standard map in accelerator 
(microtron) regime will be described. The latter model presents a unique possibility 
for quantitative study of the global critical structure. Particularly, a new part of this 
structure has been found which size was surprisingly large. Finally, in Section 5, the 
main results of the present study are summarized. In addition, the first preliminary 
empirical evidence is presented for a new regime of Poincare recurrences including 
the transition from exponential to exponential statistics. 
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2 PR in ergodic system: standard map 

Consider, first, an elementary example of ID homogeneous diffusion in momentum 
p. It can be described by a Gaussian distribution function 

where D =< (Ap) 2 > jt is the diffusion rate. Derivative fp(p,t) = df G /dp with 
boundary condition fp(0,t) = which obeys the same diffusion equation 

df _ D d 2 f 

~dt ~ ~2 ' dp^ U 
describes, then, PR to p — 0. The distribution of recurrence times (1) is simply 
related to an auxiliary function fa by 



P(r) =-Af ,„fer)* = AMO.r) = „ ^ (4) 

Here A is normalizing factor, and parameter r provides a necessary truncation of 
the preceding diverging expression at small r. It characterizes the dynamical time 
scale of the diffusion (cf., e.g., free path in molecular diffusion). If the motion in p 
is actually bounded (see below), Eq.(4) describes initial free diffusion. 

2.1 A little of theory 

Now, consider in more details another simple model - the kicked rotator - described 

by the so-called standard map: 

p = p + K sin a; (mod L) , , 

x = x + p — p 

on a torus (0<x<27r, < p < L = 2irn, n = 1,2, ...). 

We seek a solution f(p,t) of diffusion equation (3) with the boudary condition 

/(0,f) = (6) 



which provides a loss of probability because of PR to p = (and to p = L). The 
orthogonal and normalized eigenfunctions of the diffusion equation for this problem 
have the form {k > 1 is integer) 

9kip) = 1/7 sin j (7) 



with the corresponding eigenvalues 



L \ L 
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which describe the decay rate of the eigenmodes (7). In Eq.(8) the diffusion rate is 

D(K) = ■ C(K) (9) 



with the dynamical correlation function |TU 



C(K) « 1 - 2J 2 (K) + 2Jl{K) (10) 

where J 2 {K) is the Bessel function. 

The set of eigenfunctions (7) and eigenvalues (8) provides a general solution of 
the diffusion equation with boundary condition (6) for an arbitrary initial distribu- 
tion fo(p) = f{p, 0). Peculiarity of PR statistics P(r) is just in a very particular 
initial condition. Specifically, for a single trajectory in numerical experiments the 
recurrence time r is determined by the two successive crossings of the exit line which 
is, in the model under consideration, p = mod L. Hence, the initial distribution is 
concentrated right here: \p\ < K. The condition for a trajectory with initial p > 
to cross the exit line reads: p + K ■ sinx < 0. Whence, the probability of crossing is 
proportional to arccos (p/K), and the normalized initial distribution can be taken 
in the form: 

/o (k) = arccos (k ) ' 0< P^ K<L ( u ) 



An example of fo is shown in the insert to Fig.l. It is convenient to chose fo(p) on 
the one side of the exit line which is possible due to the symmetry of eigenf unctions 
(7). 

The difficulty with such an initial condition is in its narrow width which is always 
comparable with the dynamical scale (both are ~ K, a single kick). This violates 
the diffusion approximation for the exact integro-differential kinetic equation. A 
simple remedy is well known, for example, from the theory of neutron diffusion 



where the dynamical scale is the transport free path l n (see, e.g., Refs.flTTJ and JTJ, 
p. 689). A simple correction improving the diffusion approximation amounts to a 
relatively small shift of the boundary condition (6) from pi, = to pb = —al where 
I is the dynamical scale in our problem, and a ~ 1 is unknown numerical factor to 
be determined below from the numerical experiments. This implies an increase of 
the global scale: L = L + 2al while the initial distribution remains unchanged as 
it is obtained directly from the dynamics (5). Notice the corresponding change in 
eigenfunctions (7). 

The general solution of the diffusion problem is given by 

oo 

f(p,t) = E fk-gk{p)-e-^ 1 (12) 

k=l 

where the expansion coefficients fk are determined by the initial condition (11): 

[K dp [2 . (<IckK (p \\ (p 

fk = Jo KiL Sm {—{K +A ))- aiCCOS {K 

nr cos(s£A) ■ (1 - 3 {sk)) + sin(sAA) • H (s&) 
V 2L sk 

^ / l + ^V(l + ^ 2 )) (13) 



4V2L V Ti 

Here H is the Struve function, A = al/K ~ 1 (see below), and s = nK/L <C 1 is a 
small diffusion parameter. The latter approximate expression in (13) holds true for 



Now, the PR statistics is described by 



P(r) = J Q L f(p,r)dp = p /fce-^^/f^sm (^j dp 



2 ATT^ fk 



. 2L £ ^e-^ (14) 

71 m=l K 

with A; = 2m — 1 because only odd modes contribute to the integral. 

Asymptotically, as r — > oo, PR decay exponentially (Poisson statistics) 

P(r) -> | • fl + —)■ e^ 1T = F x • e" 71T (15) 



2 V 7T J 

with the characteristic time 

1 2 L 2 

n = - = - ■ — 16 

7x n 2 D 

which is determined by the first (most slow) mode m — k — 1, and which is of the 
order of the global diffusion time. 

The factor iq in Eq.(15) characterizes the share of asymptotic exponential decay 
which is small in the diffusive regime due to s C 1. The main, initial, decay is a 
power-law one. Again, due to small s, the sum in Eq.(14) can be approximately 
replaced by the integral over m to obtain: 



PW "^"^ (17) 



= 77T7ul7\ ■ 1 + — (18) 



where 

1QC(K) V ' 7r 

function C(K) is given by Eq.(10), and the approximate expression for f k in Eq.(13) 
is used. The latter is not applicable for r — > 0, so that the final expression in Eq.(17) 
is an approximate truncation of the preceding diverging relation (cf. Eq.(4)). 
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The power-law/exponential crossover time r cro is obtained from the comparison 
of Eqs. (15) and (17), and is given approximately by the relation: 

L 2 _ _ 7T 

S7TiJ lb 

Again, in the diffusive regime (L 2 3> -D) the intermediate power-law decay may be 
very long until the exponential asymptotics is reached. 

2.2 Numerical experiments 

An example of PR in ergodic case is shown in Fig.l. We use the standard map 
(5) on a torus of sufficiently large circumference L ^> K to provide a diffusive 
relaxation (s <C 1, for the opposite limit of ballistic relaxation s^l see Section 4 
below). How strange it may seem, the conditions for ergodicity even in such an 
apparently 'simple' model are still unknown ! However, numerical experiments (see, 
e.g., Ref.fl~2|1 indicate that, at least, for a particular value of the parameter K = 7 
the share of the regular domains, if any, is negligible ( ^ 10 -9 ) besides the two small 
islets (per map's period, see Section 4 below). Fortunately, their effect on PR is also 
negligible because they are related to the accelerator mode in which the momentum 
p quickly moves around the torus, so that a trajectory immediately crosses the exit 
line p = mod L (cf. Section 4). 

In Fig.l empirical data for a particular value of L = 50 it are shown which 
corresponds to 25 periods of map (5) in p. All the data were obtained from the 
run of a single trajectory over 10 7 iterations. Transition from a power law (straight 
dashed line) to an exponential (dashed curve) is clearly seen. 

For a quantitative comparison with the theory above (Section 2.1) we fix the 
dynamical parameter I = K^JC(K) = 12.1 where the value C(7) ~ 3 is used which 
has been obtained from a special numerical experiment. It considerably differs 
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Poincare recurrences in ergodic motion 




10° 2 5 10 1 2 5 10 2 2 5 10 3 



r 

Figure 1: Poincare recurrences in standard map (5) on a torus with exit line p = 
mod L = 507r: K = 7 ('quasiergodic' motion, see text), p = 0, a single trajectory 
of t — 10 7 iterations. Solid line represents numerical data; dashed lines show two 
asymptotics: a power law (17) (straight line), and exponential (15). Insert: the 
initial distribution (11), just prior to crossing the exit line. Logarithms here and 
below are decimal. 
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from the value C(7) « 1.78 according to approximate relation (10) just because of 
accelerator islets mentioned above. Since our model is a map, the minimal empirical 
recurrence time is r em = 1 instead of T t h = in a continuous theory (for example, in 
numerical data P(l) = 1). The corresponding corrections are negligible except the 
initial dependence for r ~ 1 (see below). 

Numerical data in Fig.l were fitted to Eq.(15) in the interval r = 500 — 1000 
iterations, and the empirical values of the characteristic time T\ = 125, and of 
the factor F 1 = 0.26 were obtained. The corresponding values of the correction 
parameter are a T = 2.3, and o:f = 0.71. The difference in these two values of a 
characterizes the accuracy of the correction which is rather poor because of a very 
narrow initial distribution (see Eq.(ll), and discussion around). Without correction 
(a = 0) the theoretical values would be: T\ = 69, and Fi = 0.07 which both are 
substantially underestimated. 

For a more systematic study the similar numerical data were computed for a 
number of L values specified by the integer n = L /2n. The results are shown in 
Fig.2. 

Dependence Ti(n) is well described by the uncorrected relation (16) for large n 
as expected. In intermediate region (n ~ 3) the agreement is father improved by 
the correction which provides a smooth transition to the ballistic limit (see Eqs.(28) 
and (29) in Section 4). In other words, the correction is not very important for the 
asymptotic decay rate because it is determined by the first eigenfunction which is 
only slightly disturbed, for large n, by the shift of the boundary. This is no longer the 
case for the amplitude F 1 which strongly depends just on the distorted region near 
the boundary p = 0. As a result the correction is most important for large n. The 
dependence F\{n) in the intermediate region remains unclear. For n^lO, s^l/3 
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Asymptotic exponential in ergodic diffusive motion 



PL, 




n 



Figure 2: Asymptotic exponential in ergodic diffusive case: numerical values of 
characteristic decay time T\ (16) (circles), and of the factor F 1 in Eq.(15) (crosses) 
vs. the torus size n = L /2ti. Uncorrected dependence is shown by solid straight 
lines which are transformed into dashed lines by the correction with the same average 
< a >= 1.22. Two dotted lines represent the theoretical dependence (without any 
corrections) for the opposite, ballistic, limit, Eqs. (28,29) 
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both relations, Eqs.(15) and (16), are in a reasonable agreement with the numerical 
data for the same average value of the correction parameter < a >= 1.22. 

Coming back to Fig.l, we see that the initial power-law decay is well described 
by a simple relation (17) with tq = 1 which is shown by the dashed straight line, 
and which would correspond to a = 0.66 ~ ap. 

3 PR with a chaos border: separatrix map 

Now we consider an opposite limit of essentially nonergodic system with a large 

chaos border and the critical structure. As an example we take the separatrix map 

which was studied in many papers (see, e.g., Refs. |2|, ||, §), and for which a new 

regime of PR has been recently observed |J. The latter was the main motivation 

for the present studies. We take the separatrix map in the form ||: 

p = p + sinx , 2Q . 
x = x — A ■ In (\p\) — po 

Here the motion is always strictly confined to the so-called chaotic layer: \p\ < Pb(x). 

Previously, the most studied case corresponded to big parameter A ^> 1. In this limit 

Pb ~ A, so that the width of the layer (2A) is much larger than the dynamical scale 

of the diffusion (a single 'kick') which, for map (20), is unity (cf. Eq.(5)). Besides 

the critical structure along the two borders, the average diffusion rate within the 

layer is nearly constant (see Eq.(10)): 

<D> „«VM)>„i (21) 

2 2 y J 

Hence, the initial decay of PR is a simple power law (1) which was observed, indeed, 
from the beginning [§, [|| (Section 1). The crossover time to a diferent law is given 
by a simple diffusion estimate: 

Tcro ~ ~ A 2 (22) 
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Unlike the ergodic case, the asymptotics of PR in the presence of chaos border is 
also a power law but with a different exponent v ~ 1.5. This is explained by a very 
specific critical structure near the border where the diffusion rate rapidly drops. As 
a result no trajectory can ever reach the exact border, even though it is approaching, 
from time to time, the border arbitrarily close (see Refs. [8], [L3| for details). 

An example of this well known behavior is shown in Fig. 3 (upper solid curve). A 
transition between the two different power laws (dashed straight lines) at r ~ 100 
is clearly seen in agreement with estimate (22). There is no sign of any exponential 
decay. Now, how does it appear in a similar model §|? 

The first observation is that in application to celestial mechanics (dynamics 
of asteroids) the parameter A of map (20) is typically rather small: A ~ 1 @. 
This drastically changes the structure of the layer. First of all, the layer width is 
reduced down to the size of a single kick. An example is shown in Fig.4. Hence, 
the diffusion approximation becomes inapplicable. Instead, the so-called ballistic 
relaxation comes into play which is much quicker. In other words, a slow diffusive 
motion from the exit line to a critical structure is replaced now by rapid jumps 
of a trajectory over the whole layer with some probability to get into the critical 
structure. Since those jumps are very irregular in a chaotic layer the PR are expected 
to decay exponentially. This is the case indeed as an example in Fig. 3 demonstrates 
(lower solid curve, A = 1 ). The exponential decay can be intermediate only as the 
trajectory is evetually captured into the critical structure, and the decay turns to a 
power law. Generally, the initial part of the power law is an approximate relation in 
that its exponent is not universal, and is even varying with r. In the latter example 
v « 1.1 which is rather different from v « 1.5 for the upper curve in Fig. 3. 

Another interesting and important question is how long is the intermediate expo- 
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Poincare recurrences in separatrix map 




r 



Figure 3: Poincare recurrences in separatrix map (20): solid lines represent numer- 
ical data; straight dashed lines show the power law with the exponents v indicated 
nearby; the values of parameter A = 10 and 1 are shown at the right; p = 0. Two 
cases with the same A = 1 differ by the exit line (see text). Insert: the same in 
semi-log scale. 
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nential? For the lower curve in Fig.3 it is rather long: r cro ps 150 which corresponds 
to the PR crossover as low as P cro xs 1CT 6 ! However, under different conditions 
with the same A = 1 the exponential is much shorter: r cro 50, and P cro ~ 10 -2 . 
The difference is in the exit line as shown in Fig. 4. 

In the latter case the exit line is usual: p — 0. The critical structure is determined 
by the two big islands comparable in size with that of the whole layer. This entails 
a rapid capture of a trajectory into the critical structure, and a fast transition to 
a final power law (with the local exponent v ps 1.65). The lower curve in Fig.3 
corresponds to the same A = 1 but to a different exit line: 

p ex = cos(x) (23) 

It is chosen in such a way to cut through both stability islands and, thus, to suppress 
any sticking to their critical structure. Then, the final power law is determined by 
the critical structure at the layer borders which is apparently very narrow and cannot 
be discerned 'by eye' in Fig. 4. Nevertheless, it does exist as the asymptotic power 
law of PR in Fig.3 proves. Moreover, the latter even allows us to estimate the 
size of the critical structure: its relative area (with respect to that of the layer) 
is A cr ~ 4 x 10 -4 , or the width (Ap) cr ~ 10~ 3 (see Section 4). This exponential 
transient is well fitted by the relation similar to (15) (up to r ~ T cro ) with T\ ~ 12, 
and Fi & 0.32. Both values are in a surprisingly well agreement with the uncorrected 
theory (a — 0, L ps 4.5, see Fig. 4) which gives T\ ^ 8.2, and Fi ps 0.35. Apparently, 
this is because the diffusion parameter s ps 0.7 ~ 1 is still not large enough. 

Now, we can summarize the conditions for the transient exponential in PR for 
a nonergodic motion: (i) fast, ballistic, relaxation, and (ii) a small measure of the 
regular domains. Besides, it turns out that the exponential PR allow for, at least, 
some estimates of that measure. A more quantitative study of this interesting re- 
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Separatrix map 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



X/2/T 

Figure 4: The phase space portrait of separatrix map (20) for A = 1, po — 0. 
All points belong to a single trajectory of 10 4 iterations. The straight solid line 
is the usual exit line p = while another one is the special exit line (23) which 
cuts through the two stability domains and thus blocks the contribution of their big 
critical structure to PR. 
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lation is convenient to continue with the standard map again. This is because the 
latter has an infinite series of the special values of parameter K = K n pa 2im for 
which there are well studied islands of regular motion with a simple scaling and of 
rapidly decreasing area. 

4 PR in microtron: the standard map again 

The main advantage of this microtron model is in that it is very simple, especially for 
numerical experiments, and well studied already. Here we are interested primarily in 
the domains of regular motion which exist for an infinite series of the special values 
of parameter K = K n pa 2im where n > is any integer. Within these domains 
(islands) \p\ grows indefinitly proportional to time which is the so-called microtron 
acceleration. It was well studied since the celebrated paper due to Veksler in 1944 



see, e.g., [14] and references therein). However, in the present paper, as well as in 



Ref . ||14|| , the main object for study is not the regular acceleration itself but rather 
the chaotic motion outside the microtron islands which is generally affected by the 
critical structure at the island borders. A picture of this scale-invariant border is 
shown in Fig. 5a in dimensionless variables 

x s = (x - x ) ■ K , p s = (p - po) ■ K (24) 

where p is a parameter of map (5), and 

K-smx = 27m, K 2 = o 2 + (27m) 2 , a = K-cosx , -4 < o < (25) 

The latter inequalities determine the stability region around a fixed point 

±xo, Po mod 2tt. In Fig. 5a and below a = —2 (the center of stability). For each 

integer n there are two islets per phase space bin 2tt x 2tt one of which is presented 
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in Fig.5a. The picture shows a single trajectory of 5000 iterations. During this 
time interval the trajectory is sticking to the critical structure very close to the 
exact chaos border which results, under particular conditions (see below), in an 
asymptotic power-law decay of PR (cf. Fig. 3 above). The island relative area (with 
respect to that of the phase-space bin) is given also by a dimensionless relation [TJj : 



A n ■ K 2 n = A {<t) « 0.17 (26) 

where the latter value corresponds to a = —2. This area rapidly decreases as island's 
number n grows. Yet, for any n — > oo it determines the asymptotic PR decay, as 
we shall see below. 

In Fig. 5b another, much smaller, microtron island is shown for comparison. In 
this case an outside, and much longer, trajectory was used which cannot ever cross 
the chaos border and enter the island. Its area is given by the same estimate (26) 
with Ao(-3.1) « 0.0038. 

The main difficulty with the microtron model for our purposes here is the rapid 
growth in \p\ within and around the chaos border. This destroys any long sticking 
of a trajectory whatever the exit line for PR (cf. Section 3, Fig. 4). To overcome 
this difficulty we used the following method. First, we have chosen the exit line in 
such a way not to cut any island. It was done simply by fixing parameter p = 
7t(t^ mod 27r) in map (5) without any change in the configuration of map's torus. 
Second, we compensated acceleration by adding the term 27m to the first Eq.(5). 
This helps, of course, for one island of each pair only. 

Now, we need to provide the ballistic regime of relaxation that is a sufficiently 
large parameter s = nK/L (Section 2.1). It is convenient to take L = 2nn, so that 
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(a) 



(b) 




-5 -4 -3 -2 -1 



1 2 3 4 5 




\j — ■ ■ 

-0.5-0.4-0.3-0.2-0.1 0.0 0.1 0.2 0.3 0.4 0.5 



Figure 5: Universal border of microtron islets in scaled variables (24): (a) - in the 
center of stability interval, o = —2, single trajectory of 5000 iterations stuck to the 
chaos border; (b) - near the edge of stability, o m —3.1, revealed by an outside 
trajectory of 10 7 iterations, K = 7. 
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the parameter 



s = ^ » (27) 
27m v y 



is nearly independent of n except a few small values of the latter. 

Neglecting any dynamical correlations of the motion (particularly, those caused 
by the presence of small microtron islets including the compensation of acceleration) 
it is straightforward to calculate the probability w, per map's iteration, for a trajec- 
tory to stay within the torus without crossing the exit line. As is easily verified, it 
is given by the relation: 



w = J^rk = dx{L ~ K - Sinx) = l{ Xm ~ ^ (1 - COSXm) ) 

2 



- « 0.363 (28) 

7T 



where 

= ( . >** (29) 

{ 7T/2 , S < 7T 

These general relations were used in Section 2.2 (Fig. 2) to draw the ballistic approx- 
imation. 

The latter expression in Eq.(28) corresponds to the value s = n used in numerical 
experiments. Without additional shift Ap = 2nn discussed above the average time 
of the exponential decay would be 

< r >= ~ ► 0.988 (30) 

m w 

For s = n the shift increases w and < r > up to 

w 1 1 

w = — + - = 1 w 0.682 

2 2 7i 

< t > = --— w 2.61 (31) 

m w 

Now we can turn to numerical experiments with this microtron model. 
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4.1 PR in microtron: numerics 

The main results of numerical experiments are presented in Fig. 6, and in the Table 
below. In Fig. 6 the points show numerical data computed from a single trajectory 
(for each n) up to 3 x 10 11 iterations (for the largest n = 5000). The straight solid 
line is the fitted intermediate exponential with the decay time < f >= 2.41 in a 
good agreement with the expected theoretical value 2.61 in Eq.(31). This justifies 
neglecting dynamical correlations assumed in the above theory in ballistic regime. 

The exponential/power-law crossover time systematically increases with n that 
is with the decrease of the microtron island area (see Table). The power-law tails 
of PR were fitted by the expression 

Pn(r) = ^ (32) 

Remarkably, all values of the exponent were found to be close: v « 2. The relation 
of this expression to the size of the critical structure is based on the following 
hypothesis: dependence (32), fitted to the tail of PR, can be extrapolated back to 
r = 1. If true, it allows us to interpret the parameter A cr (n) as the relative area of 
the whole (global) critical structure around the corresponding microtron island of 

One could expect that both areas are comparable: A cr {n) ~ A n . Surprisingly, 
this is not the case (see Table, third column; the data in fourth column will be 
discussed below). Their ratio R = A cr /A n ~ 100 is not only very large but also 
slowly increasing with n according to the following approximate empirical relation 

R{n) « 50n 1/4 (33) 

The origin of this small correction to a simple scaling R m const remains unclear. 
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Poincare recurrences in microtron 




T 

Figure 6: Poincare recurrences in microtron (points) for the values of the microtron 
parameter n given at the right. Straight line shows the intermediate exponential 
with the fitted decay time < f >= 2.41 (cf. Eq.(31)). Solid curves give the fit of 
asymptotic power law (32) with nearly the same exponent v rs 2 extrapolated back 
to r = 1 (see text). 
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Table. Global critical structure around microtron islets 
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7.7 



In any event, the size of the whole critical structure seems to be much larger than 
expected. This main outer part of the structure looks ergodic, and forms a sort of 
halo around the usually narrow inner part with a typical admixture of chaotic and 
regular components of motion. The former reminds the ergodic critical structure 
around a parabolic fixed point, that is the limiting case of an island of zero size, 
studied in Ref . |T5| . In a sense, such a halo is some 'hidden' critical structure, without 
internal chaos borders but with apparently strong correlations in the motion which 
keep a trajectory within this relatively small domain. 

Now, the principal question to be aswered reads: is the observed halo a real 
physical structure or the result of a wrong interpretation of the empirical data using 
the above extension hypothesis? 
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4.2 Exit time statistics 

To clarify this question a new series of numerical experiments was undertaken. To 
this end, the exit times from the halo, instead of recurrences, were measured. Such 
a method was recently successfully used in the studies of the critical structure in 
Ref . 1 16fl . In the problem under consideration here the measurement of exit times 



was organized as follows. A number (typically 100) of trajectories with the initial 
conditions homogeneously distributed over the circle around a microtron island (see 
Fig. 5a) were run until they leave the interval (0 < x < it). The dependence of 
the average exit time r ex for a series of the circles with increasing radius p s as a 
function of the area within a circle A s = 7rp 2 s (in scaled variables (24)) was thus 
computed. The minimal circle of radius p s = 3 touches the island, and comprises 
the area A m i n = 28.3 while island's area in these units is A sn = 6.72, the minimal 
ratio being R ex = A s /A sn = 4.21. 

The main results of this measurement are shown in Fig. 7 for 8 different values of 
parameter n up to n = 10 7 with the island small as A n « 4 x 10" 17 ! This 

is completely out of reach for the PR method (cf. Ref. ||16||). The difference is in a 
rather short exit time from the halo, we are interested in, as compared to the long 
recurrence time on the tail where it is eventually separated from the exponential 
(Fig.6). 

The main result revealed in Fig. 7 is a transition between the two different scal- 
ings. One, for relatively large r ex , is the standard critical scaling shown by the 
dashed line which is the fitting of numerical data to the relation 

R ex (r ex ) « ^ill (34) 

with R ex (l) ~ 50. As expected, this part of the data does not depend on n. More- 
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Scaling of exit times from halo 
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Figure 7: Scaling of exit time from halo: the ratio R ex of the critical to island area 
vs. the average exit time out of the former (points) from a hundred trajectories (for 
each of 8 values of n) with initial conditions on a series of circles (see Fig.5a and 
the text). Dashed line shows the standard critical scaling (34) which breaks down 
at some values of crossover time r cro depending on n (insert). 
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over, scaling (34) is in a good agreement with the PR tail in Fig. 6. The relation 



between the two is well known [g, ||, [L4], [1d[| . Generally, the power-law PR statistics 
is descibed by (cf. Eq.(32)): 

W * m „ Mt) <i> ( 35) 

T u T 

where < r > is the average PR time, and the latter expression is obtained from the 
ergodicity within the chaotic component of the motion. Using approximate relation 
t 2 r ex |L6| we obtain 



where i? = A cr /A n ps A(l)/A„ (see above). For integer map's time 

<T>= I T (^-F^) =<M (37) 

where £ is the Riemann function. Whence, for z/ = 2 the relation (36) gives (see 
Eq.(34)) 

i? « —R ex {l) « 160 (38) 

which is in a reasonable agreement with numerical data (third column in Table). 
However, unlike the data in Fig. 6 where the actual power-law scaling is not seen 
under much larger exponential transient, the data in Fig. 7 clearly demonstrate that 
the critical scaling does not reach the limit r = 1 assumed above. Moreover, the 
crossover time r cro increases, and hence the size of the global critical structure (R ex ) 
decreases, as n grows (Fig. 7, insert). The increase of r cro must have an upper 
bound because otherwise the critical structure near the chaos border would be also 
destroyed in contradiction to the detailed studies of that in anomalous diffusion ji~4 



Indeed, the empirical dependence r cro in Fig. 7 (insert) can be fitted reasonably well 
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by the expression 



4 66 

6-69 - ^ (39) 



The upper limit in r cro corresponds, according to Eq.(34), to the lower limit in the 
ratio R ex > 7.47. Combining Eqs. (34) and (39) we obtain approximately 

R ex (n) « ( 4 °) 

The empirical values of this dependence are given in Table (fourth column). They 
indicate much smaller, yet still a fairly large, size of the critical structure as compared 
to the limiting estimate for PR (third column). The former seems to be more reliable 
and realistic. A different, new and unknown, scaling in Fig. 7 for r ex < r cro requires 
further studies. What is of importance here is the termination of the critical scaling 
at a finite r ex = r cro . This determines the outer border of the critical structure. 

5 Discussion: a new puzzle 

The original motivation of these studies was the unusual exponential transient ob- 
served in PR in the presence of chaos border ||. However, in the course of in- 
vestigating the mechanism and conditions of this phenomenon a more interesting 
observation has come out. It suggests the existence of a new, unknown to my knowl- 
edge, part of the critical structure surrounding, like a halo, the well-known inner 
part close to the chaos border. In spite of some contradictory empirical evidence 
the halo apparently occupies the most of the global critical structure. In any event, 
in the microtron model considered in this paper the area of the halo is much larger 
than that of the regular island inside it, even according to the minimal estimates 
(see Table and Fig. 7). 

As is well known, the scaling of the peripheral part of the critical structure is 

27 



generally nonuniversal, at least quantitatively, in the sense of the corresponding 
power-law exponents, for example However, it might be nevertheless typical 
qualitatively as it appears in our model. In this respect, it would be interesting to 
look at different examples of the global crtical structure. One possibility is to use 
the same model with a fixed parameter L = 27m (Section 4) but for different values 
of the stability parameter a in Eq.(25). First preliminary numerical experiments 
have been done for 9 values of a within the whole stability interval (—4 < o < 
0; 2tx < K < 7.45) including the 'quasiergodic' case K = 7 used in Section 2 for 
other purposes. In all cases but the latter the PR behavior was similar to that in 
the main series of numerical experiments (Fig.6), at least qualitatively. However, 
just for K = 7 a sudden surprise has emerged which is presented in Fig. 8. 

In spite of a very long run (10 11 iterations) no clear sign of the expected power- 
law decay is seen. A small deviation from the final exponential at the end of the 
dependence is a typical feature due to a poor statistics (cf., e.g., Fig.6). The first 
exponential is close to the expected one with the fitted decay time < f >= 2.41 as 
compared to the theoretical < f >= 2.38 (see Section 4). For the second exponential 
the empirical decay time < f >= 23.1 is about 10 times longer. This means that a 
trajectory is kept within (sticks to?) a certain domain but not in a way it does so in 
the usual critical structure. Moreover, the relative area ~ 2 x 10~ 4 of this peculiar 
domain, estimated similarly to A cr (n) in Eq.(32), is small and is comparable with 
that of the island inside (Fig.5b): A 7 w 7.8 x 10~ 5 . This island does have a chaos 
border, yet contrary to usual behavior, it does not produce any appreciable power- 
law decay of PR. Another preliminary remark is that a more careful inspection 
of Fig. 5b seems to suggest a different, more regular than usual, structure of the 
chaos border for K = 7 (cf. Fig. 5a). Certainly, this 'anomaly' deserves further 
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K7 puzzle 
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Figure 8: Poincare recurrences for L = 2ir, and K = 7. The solid line shows 
numerical data from a single trajectory of 10 11 iterations. Two dashed lines are 
fitted exponentials with the average decay time < f >= 2.41, and < f >= 23.1, 
respectively. 
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investigation. 
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